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We congratulate the authors for their overview 
paper discussing modeling techniques for spatial ex- 
tremes. There is great interest in spatial extreme 
data in the atmospheric science community, as the 
data is inherently spatial and it is recognized that 
extreme weather events often have the largest eco- 
nomic and human impacts. In order to adequately 
assess the risk of potential future extreme events, 
there is a need to know how the characteristics of 
phenomena such as precipitation or temperature 
could be altered due to climate change. 

Because of the high interest level in the atmo- 
spheric science and (more broadly) the geoscience 
communities, it is imperative for the statistics com- 
munity to develop methodologies which appropri- 
ately answer the questions associated with spatial 
extreme data. Davison, Padoan and Ribatet (2012) 
provide a comprehensive overview of existing tech- 
niques that can serve as a useful starting point for 
statisticians entering the field. That the paper is 
written as a case study helps to illustrate the advan- 
tages and disadvantages of the various methods. We 
hope that this Swiss rainfall data will serve as a test 
set by which future methodologies can be evaluated. 

The authors analyze data which are annual max- 
ima. This is natural from the classical extreme value 
theory point of view whose fundamental result es- 
tablishes the limiting distribution of Y = ( VT=i 
■ • • j VT=i -^Di) T to be in the family of the multivari- 
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ate max-stable distributions. In practice, modeling 
vectors of annual maxima seems less than ideal, and 
it is not clear how much dependence information 
is lost by discarding the coincident data. Scientists 
in other disciplines can be uncomfortable with the 
idea of constructing data vectors of events which 
most often occur on different days. We are aware 
that there is current work to extend spatial extremes 
work to deal with threshold exceedances, and we 
look forward to that work appearing in the litera- 
ture. 

Davison, Padoan and Ribatet (2012) divide the 
spatial approaches into three categories: latent vari- 
able models, copulas, and max-stable process mod- 
els. In Section 7 they do a very nice job of detail- 
ing the strengths and weaknesses of the three ap- 
proaches. However, it seems that the article does 
not make clear enough that the aim of the latent 
variable approach is fundamentally different than 
the aim of a copula or max-stable process model. 
As the authors state in Sections 2.2 and 2.3, cur- 
rent modeling of multivariate (or spatial) extremes 
requires two tasks: (1) the marginals must be esti- 
mated and transformed to something standard (e.g., 
unit Frechet) so that (2) the tail dependence in the 
data can be modeled. The latent variable model is 
a method for characterizing how the marginal dis- 
tribution varies over space, that is, task 1. In con- 
trast, both copula models and existing max-stable 
process models explicitly model the tail dependence 
in the data once the marginals are known, that is, 
task 2. We refer to the dependence remaining after 
the marginals have been accounted for as "residual 
dependence," as Sang and Gelfand (2010) described 
the random variables after marginal transformation 
as "standardized residuals." 

Davison, Padoan and Ribatet (2012) are correct 
to point out (Figure 4) that using a latent vari- 
able model is inappropriate for applications where 
the joint behavior of the random vector is required. 
However, there are applications which aim only to 
model the marginal behavior. There is a long history 
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of producing return level maps such as those shown 
in Figure 3 of the manuscript. For instance, the re- 
cent effort to update the precipitation frequency at- 
lases for the US (Bonnin et al., 2004a, 2004b) aimed 
only to characterize the marginal distribution's tail 
over the study region. Bonnin et al. (2004a, 2004b) 
employed a regional frequency analysis (Dalrymple 
(1960); Hosking and Wallis (1997)), approach which, 
like the latent variable model approach, aims to bor- 
row strength across sites when estimating marginal 
parameters. As Davison, Padoan and Ribatet (2012) 
clearly show, explicitly modeling residual depend- 
ence requires considerable effort, and when only the 
marginal effects need to be described, we feel it can 
be appropriate to ignore the residual dependence so 
long as one recognizes the limited scope of the ques- 
tions that such an analysis can answer. 

In situations where the joint behavior of multi- 
ple locations must be described, then one must ex- 
plicitly model the residual dependence. As Davi- 
son, Padoan and Ribatet (2012) show, dependence 
models not specifically designed for extremes may 
be inadequate to capture tail dependence. However, 
models such as the extremal copulas or max-stable 
processes do not easily lend themselves to current at- 
mospheric science applications with hundreds, thou- 
sands, or tens-of-thousands of locations. There are 
obvious avenues to explore toward adapting the pair- 
wise likelihood methods for large spatial data sets, 
but, to date, pairwise likelihood methods have only 
been applied to applications similar to the one in 
Davison, Padoan and Ribatet (2012) with roughly 
50 locations. We imagine scaling the methods to the 
size of current applications will be nontrivial, and 
perhaps new inference procedures or more computa- 
tionally-feasible extremes dependence models will 
need to be developed. Until appropriate extremes 
techniques are available, people will continue to be 
tempted to apply high-dimensional models devel- 
oped to describe nonextreme data (e.g., a Gaussian 
copula) to model tail dependence. 

Most of the spatial extremes work to date has 
been primarily descriptive in nature. Such analy- 
ses are useful in assessing risk (i.e., the probability 
of an extreme event), but do not help to explain 
the underlying causes of extreme events. There is 
a desire in the atmospheric sciences to move beyond 
descriptive analyses and toward analyses which en- 
hance understanding of the processes which lead to 
extreme events. For example, Sillmann et al. (2011) 
establish a link between extreme cold temperatures 
in Europe and a blocking phenomenon in the North 



Atlantic, Maraun, Osborn and Rust (2011) link ex- 
treme precipitation in Europe to large-scale airflow 
covariates, and Weller, Cooley and Sain (2012) link 
extreme precipitation on the Pacific coast of North 
America to surface pressure patterns. Since it is gen- 
erally believed that climate models are better at rep- 
resenting processes at large-scales, establishing links 
between extreme events and large-scale phenomena 
enable one to better conjecture how the nature of 
extreme events will change with the climate. While 
none of the analyses cited above involved extensive 
spatial modeling of extremes, it is foreseeable that 
science will move in this direction. 

Finally, undertaking a pairwise likelihood fitting 
of a max-stable process model is challenging and 
would be beyond the capabilities of most geoscien- 
tists. The authors are to be commended for develop- 
ing the SpatialExtremes (Ribatet (2011)) package 
in R which enables the general scientific community 
to utilize these methods. 
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